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Abstract 


The dynamic characteristics of rotor blades play a significant role in the 
overall performance and stability of helicopter. A time domain integration, 
approach is used to solve the set of non-linear differential equation defining 
behavior of rotor blade / fuselage with non-dimensional time. 

Using Multiblade Coordinate Transformation an attempt to transforms the 
rotating coordinates into a non-rotating frame is made in order to study the 
time response of rotor degrees of freedom. 

An attempt is made to develop a parameter, which measures the rate of 
divergence (instability) or convergence (stable response) and its level of 
dependence on the operating and/or initial condition. 

Two sets of table are also drawn in order to study the effect of support 
damping and with zero support damping on the stability (or instability) of 
the system. 
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CHAPTER 1 


INTRODUCTION 


1.1 Bearingless Helicopter Rotor Blade: 

The main rotor system of a helicopter is the most important component of the vehicle. 
The rotor system has to fulfill multiple roles of generating lift, propulsive and control 
forces. Therefore, the dynamic characteristics of rotor blades play a significant role in 
the overall performance and stability of the vehicle. There has been a continued effort 
to develop a rotor blade and hub configuration which is mechanically simple yet 
efficient. 

The development of a bearingless rotor system aims to eliminate the 
conventional pitch bearing as well as the external damper by incorporating a 
specialised elastomer with high loss factor. In this rotor system, the blade is attached 
to the hub through a flexible structural element called flexbeam. The flexbeam is 
designed to provide the required stiffness in the flap and lag bending deformations of 
the blade, but it is highly flexible in torsion. Surrounding the flexbeam, there is a stiff 
cuff denoted as torque tube which is attached to the blade-flexbeam junction at the 
outboard end and to a pitch link at the inboard end as shown in figure [1.1(a)]. The 
pitch control of the blade is achieved by rotating the torque tube through up/down 
movement of the point P, which in turn twists the flexbeam. An elastomer is placed 
between the torque tube and the flexbeam to provide adequate lag damping. It also 
serves the purpose of a spacer between the torque tube and the flexbeam. Though a 

f 

bearingless rotor is mechanically simple , see figure [1.1(b)], its dynamic analysis 
becomes very complicated due to the presence of multiple load paths and non-linear 
elastomeric damper. It is important to point out that in a practical bearingless rotor 
blade configuration, the flexbeam and blade properties will be different. In this thesis, 
it is assumed that the blade arid the flexbeam have same uniform properties. The 
reason for this assumption is (i) practical blade data is not available (ii) by considering 



uniform property for bearingless rotor, one can obtain good physical insight of the 
effects of non-linear elastomer. 

1.2 Background Information: 

Over the years, considerable effort has been put in by the helicopter industries in 
developing soft-in-plane bearingless main rotors with acceptable performance and 
stability characteristics. Huber[l] presented a comprehensive review of the 
bearingless main rotor programme pursued at different helicopter industries in the 
world. A brief history on the development of bearingless rotors is provided in Refs. 
[2,3]. Lockheed developed the first prototype of a bearingless rotor in mid 1960’s 
having “matched stiffness” (equal stiffness in both in-plane and out of plane bending). 
During test flight, the prototype encountered an air-resonance instability which was 
eliminated by increasing the fuselage inertia. During late 1970’s intense activity 
continued in developing the prototype of Bearingless Main Rotor (BMR). The Boeing 
Helicopter company developed a BMR with a high fundamental flap frequency and 
tested it on and MBB BO- 105 Helicopter. It was observed that the landing gear of B0- 
105 had to be stiffned in order to improve the ground resonance stability. Another 
design was proposed by the Bell Helicopter with a four bladed M680 rotor. An 
elastomeric damper was provided to augment edgewise damping. In 1990 the Boeing- 
Sikorsky company developed a five bladed BMR (for the new helicopter known as 
RAH-66 Comanchi) with elastomeric dampers. This vehicle is currently under going 
tests, [2]. In Europe, Eurocopter recently built and flight tested EC135 helicopter with 
a BMR, [4]. For the analysis and design of bearingless rotor blades, most of the 
rotorcraft manufacturers have developed their own in-house methodologies- and 
computer codes. Hence, very little information is available about them in the open 
literature. Recently Pohit,G. et.al [5] have analyzed the influence of non-linear 
elastomer on isolated lag dynamics and rotor/fuselage aeromechanical stability and 
studied the non-linear free vibrations of the rotating blade under lag bending. 
Numerical perturbation technique is used to determine the frequency-amplitude 
relationship. A conclusion has been drawn that up to a fairly high value of amplitude, 
a seventh order expansion for the resorting force is sufficient to correctly predict the 
frequency-amplitude relationship. They have also carried out amplitude dependent 
stability analysis for a coupled rotor/fuselage system under ground resonance 



resonance condition. The non linearities of elastomer were treated by an equivalent 
linearisation approach. They observed that the elastomer amplitude does not alter the 
stability of the system in the region of maximum instabilities, but the stability of the 
progressive lag mode increases appreciably with increasing amplitude. Also, the 
location of the elastomer and torque tube have significant influence on the stability of 
the system, there is an optimum location for the torque tube providing least instability 
in the system. 


1.3 Elastomeric Damper: 

Elastomeric materials are becoming increasingly popular for vibration attenuation 
applications due to the tremendous advantage they offer over mechanical dampers. In 
the case of bearingless rotors, conventional hydraulic dampers are being replaced by 
specialised elastomeric material with high loss factor. Under dynamic conditions, 
elastomeric materials exhibit visco-elastic behaviour dissipating energy through 
hysteresis. The characteristics of an elastomeric material are also non-linear with 
respect to the amplitude of motion. Consequently it has been a challenging task to 
develop an analytical model that can accurately represent the behaviour of the 
elastomeric dampers. 

The non-linear characteristics of elastomers have been analysed in detail in 
several recent publications. In Ref. [5], the non-linear characteristics of the 
elastomeric material have been idealised by a combination of linear and non-linear 
restoring and dissipative elements. Figure (1.2) represents an idealized model 
consisting of a non-linear spring, coulomb damper and hysteric damper, which differs 
from the model proposed by Gandhi et. al. [7] in the sense that their model has a 
viscous damping element. Figure (1.3) shows the variation of in-phase (k') and 
quadrature (k") stiffnesses, which correlate well with the experimental data. 

The model developed in Ref. [5] has been used in the present study. 


1.4 Aeromechanical Stability Analysis of Helicopter: 

The aeromechanical stability of a helicopter is a complex phenomenon, involving 
coupling between the rotor and body degrees of freedom. The rotor lead-lag 



regressing mode usually couples with the body pitch and roll modes to cause 
instability. This phenomenon is denoted as ground resonance when the helicopter is 
on ground. The aeromechanical stability is severe in nature. In case of bearingless 
rotors, a specialised elastomer with high loss factor is used as a damper. 

The analysis of aeromechanical stability characteristics of a coupled 
rotor/fusclage system requires the development of a mathematical model representing 
the dynamics of the coupled system. Over the years, considerable progress has been 
made in modelling the system dynamics. Friedmenn and Venkatesan [7,8] developed 
analytical models for the aeromechanical stability of a coupled rotor/fuselage system 
and compared their theoretical results with the experimental data presented 
Venkatesan [10] analysed the effects of pitch-lag, pitch-flap and structural flap-lag 
couplings individually and in combination, on the stability of lag mode under ground 
resonance condition. In Ref. [6], ground resonance stability, using an idealised model 
of the coupled rotor/fuselage system has been analysed and is shown in figure [1.4]. 
In this model, the blade is assumed to undergo only lag motion and the elastomer is 
represented by non-linear spring and damping elements. The torque tube is assumed 
to be rigid and massless. The fuselage is allowed to undergo rigid body pitch and roll 
motions. Emphasis is placed on the effect of non-linearity of the elastomer on the 
ground resonance problem. The nonlinearities were linearised for a given amplitude 
using equivalent linearisation approach performing a linearised stability analysis. He 
observed that effect of amplitude on ground resonance analysis is that the progressive 
lag mode becomes more stable with increasing amplitude for the baseline 
configuration (^=0.1 and <^ 2 =0.25) and C x =0, C y =0. 

The present study is an extension of the study carried out in Ref. [6]. 


1.5 Objective of the Present Study: 

The major objectives of the present study are to: 

(i) Obtain equation of motion for coupled rotor/fuselage system including 
elastomer nonlinearities, under ground resonance condition. 

(ii) Analyse the amplitude dependent aeromechanical stability by directly 
integrating the governing differential equations in time domain, for various 
sets of operating and initial conditions. 

(iii) Identify the region of instability. 



(iv) 


Formulate a measure of instability of the system and study its variation with 
both operating and initial conditions. 





(b) Actual rotor hub (taken from Ref. {4]). 

Fig 1.1 Bearingless rotor hub and Blade configuration 
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Fig. 1.2 Idealised elastomer model 
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Fig. 1.4 Coupled rotor / fuselage model 





CHAPTER 2 


EQUATION OF MOTION 


2.1 Ground Resonance Analysis 

The aeromechanical stability of a helicopter is a complex phenomenon, involving 
coupling between the rotor and the body degrees of freedom. The rotor lead-lag 
regressive mode usually couples with the body pitch and/or body roll to cause the 
instability. When the helicopter is on the ground, such instability is referred to as 
ground resonance. In the case of bearingless rotor, the ground resonance is usually 
avoided by providing elastomeric damper to the rotor blade and damping in the 
landing gear. 

The formulation of equation of motion presented in this chapter is for coupled 
rotor/fuselage system under ground resonance condition. 


2.2 Bearingless Rotor Model: 

An idealised model of the coupled rotor/fuselage system is shown in figure [1.4]. The 
blade is assumed to undergo only the lag motion. The elastomer is represented by 
non-linear spring and damping elements. The torque tube is assumed to be rigid and 
massless. The effect of pitch link in the lag mode is neglected. The hub is located at a 

height ( h ) above the C.G. of the fuselage. The fuselage is allowed to undergo rigid 
body pitch and roll motions. The characteristics of the fuselage support (landing gear) 
are represented by a linear spring and a linear damper. 

2.3 Formulation of the Elastomeric Damper Model: 

The non-linear characteristic of an elastomeric material can be idealised by a 
combination of linear restoring and dissipative elements, since the static stress-strain 



curve of an elastomeric material is non-linear. Figure [1.2] represents an idealized 
model consisting of a non-linear spring, a coulomb damper and hysteric damper. The 
reason for choosing a coulomb damper is that the damping force is very high at low 
amplitude. In addition, coulomb and hysteric dampers provide damping forces which 
are independent of frequency. This type of idealization differs form the model 
proposed by Gandhi et. al [30], in the sense that their model has viscous damping 
element. 

The non-linear force-deformation relation of the spring is assumed to be of the 

form 

Fj = K \X — IC 3 x 3 K §x^ — K.’jx'i (2. 1 ) 

Where F s is the force exerted by the restoring element under deformation x with Ki, 
K 3 , K 5 , and K 7 as positive constants. The constitutive differential equation of the 
elastomer model under external loading can be written as 

K\X- K 3 x 3 + K s x 5 - K 1 x 1 + F sgn | x | +— x = Externally applied force 

£0 

Numerical values of the system parameters of the elastomeric damper model are taken 
directly form Ref. [5] and is presented in table (2.1) 

2.4 Formulation of Equation of Motion 

Several assumptions have been made while formulating the equation of motion. These 
are: 

i. The rotor blade is an Euler-Bemuolli beam made of homogenous isotropic 
material, and undergoing lag deformation. 

ii. The torque tube and pitch link are rigid and massless elements (usually, the 
torque tube has a very high stiffness relative to that of the flexbeam or rotor 
blade section). 

iii. Aerodynamic effects are not considered since emphasis is placed on the effect of 
non linearity of the elastomer on the ground resonance phenomenon. 

iv. Effect of pitch link is neglected. 

The dynamical equations of motion of the coupled rotor/ fuselage model consist of 
two sets of equations : one corresponding to the blade motion and the other 



representing the motion of the fuselage. The equations of motion are derived using 
Hamilton'’ principle. 


2.4.1 Blade Equation 

The important ingredients required for the formulation of blade equations are: 

(a) The kinetic energy of the blade. 

(b) The strain energy of the blade. 

(c) The strain energy associated with the elastomer. 

(d) Dissipative energy associated with the elastomer. 


2.4.1.1 Kinetic Energy of the Blade: 

The kinetic energy of the k th blade in an N-bladed rotor system can be written as 


no= 


du _ _ dx dy . 

Qv + — cosy* +— smy* 

i dt dt dt 


dv 


_ dx . ay I 

+ Qx+Qu smy* + — cosy* V 

dt dt dt 


dy 


\ 2 


J 


( 2 . 2 ) 


dx 


Where x(~Q y h ) and y(=-6 x h ) are the perturbational displacements of the hub 

centre due to fuselage motion in pitch and roll models, respectively. \]/k is the azimuth 
angle of the K blade, u and v are the axial and in-plane displacements of the blade. 


2.4.1.2 Strain Energy of the Blade due to Bending: 

The Expression for strain energy due to bending is 



2.4.1. 3 Strain Energy Associated with the Elastomer: 

The expression for the strain energy of the spring due to deformation A/ is 


U, 


where A/ = 


v(L 2 ,i)-(L 2 -A) gV(( t 2 ' <) -v(A.O 

OX 


(2.4) 



2.4.1. 4 Dissipative Energy Associated with the Elastomer: 

The expression for the virtual work done by the dissipative (coulomb and hysteretic) 
elements of the elastomeric damper in the lag mode is given by 

5W d =F d 5 A, (2.5) 

where A/ is the deformation of spring in lag mode 

= v(l 2 , <) - (l 2 - L , ) - vtt, , 0 

ox 

( 2 . 6 ) 

and 

F d =-^A l +Fsgn\A,\ 


(2.7) 

where h is the hysteretic damping coefficient 
F is the coulomb damping coefficient 

03 is the frequency, experimentally data fitted value, value taken is 3.3Hz. 


2.4.2 Fuselage Equations 

The in-plane (S v ) and the radial (S r ) forces at the root of the k th blade are given by 

l L 

S v = j (-mv + mQ 2 v)dx-(xsin x i , k +ycos x ¥ k ) Jmdx (2.8) 

0 0 

and 


L L L 

S r = ^mQ. 2 xdx- j2Qv/m£c-(;ccos^£ +ysin\i/j c )jmdx (2.9) 

0 0 '0 

In-plane (S v ) and radial forces (S r ) can be re-written in non-dimensional form as 
given below; 


In-plane force: 

l _ 

• S„ = -mn 2 i 2 [ g )i| ] \y, t (t) ■ - (x sin - ? cos 4-* ) • (mi) 

0 



Since 


x = hQ y and y = -h & x 

i _ 

S v =-mQ 2 L 2 [j(fn^)d^]y lk (z)-h -(mL)-{6 y sinT* +6 X cosT*} (2.8.a) 
0 


Radial force: 

1 1 - 

S r =-(n,n 2 L 2 )-2(mCl 2 L 2 )[ \^)d^]y Xk {x)-h -{mL)-{Q y cos** -0* sinT'*} 

2 o 

(2.9.a) 

The forces acting on the hub can be obtained by resolving the blade- root shear forces 
in the non-rotating frame and summing over all the N-blades, in our case N=4. 

In X-direction: 

PH , = 2 ;(cos'P*S r -sin ¥*$„) 
a:=i 

On substituting the values from eq 2.8 (a) and 2.9(a) we get 

l _ 

P^,=(- 2 «Q 2 L 2 ).[J^ 1 (^].[C 0 S'F, ^,+ 008^2^2+008^3^3+008^4^4 

0 

1 _ 

+ (mO 2 ! 2 ) • [ (g)d$] ■ [sin'F, y x { + sinTj y 12 + sin ¥ 3 y n + sin ^4 y u ] 

0 

In Y-direction: 

4 

p H y = J](sin'T A 5 r - cos V k S v ) 

K = 1 

On substituting the values from eq 2.8(a) and 2.9(a) we get 

1 _ 

PH, =(-2mQ J i 2 )-[^ 1 (|)d|].[sm>P 1 7 M +sin>P 2 ^j+sin'i'j y 13 +sinT 4 y lt ] 

0 

1 ___ 

-(mQ 2 X 2 )-[J^ 1 (^)^]-[cos v P 1 y 11 +cosT 2 y 12 +cosV 3 y l3 +cos'F 4 y 14 ] 
0 

Considering the support structure as a combination of mass-spring-damper system 
having pitch and roll motions, equations of motion of the support can be written on 
balancing the moments in pitch and roll direction as 



I9 X + C6 X +K9 X --PHy -h 


and 


I9y+C9y+K9 y - PH x -h 
On substituting the expressions for PH x and PH y the 
can be written as 

I 

Jij 

X 0 X +K X 0 x = 2 (— ){sin v F, ^ n + sin^ 2 sin V F 3 y 13 + sin v F 4 y l4 } 

hJ . 


K9 V - PH x • h 

fuselage equation of motion 


I,0 X +C 

+ (— ){c°S v Fi y l{ -l-COS^ > ; 12 + cos '^3 ^13 + 
L 


and 

Iy9y+Cy9y+Ky9y 


cos^ y 14 } 

( 2 . 10 ) 


+ cos^ y l2 + cos ¥3 j / 13 + cos ¥4 y 14 } 
+ (y ) {sin ^ y ii + sin T 2 y 12 + s 


K V 9 V =-2(y){cos'F 1 >; 11 


-sin^ j> 13 + 


sin^ y x 

(2.1 


where, 

I x , I y are non-dimensional mass XXXUXXXVXXIO ' 
l and pitch modes, respectively. 


moments of inertia of the support structure 


the C.G. in roll 

I / = — - — 

*' y ml?h 

C x , C y are non-dimensional linear damping coefficients of the support 
roll and pitch mode, respectively. 

C C = — 

*’ y mQL 2 h 

K x ,K y are non-dimensional stiffiiesses of the support structure in roll 
modes, respectively. 

and 


1 _ 

o 4 =J=|A(«^ 


0 



4* is the azimuth location of the K-th blade. Locations of all the four blades, for 
rotor/fuselage model, are given by 

V k = V + K=l,2,3,4 (2.12) 


first blade 
second blade 
third blade 
fourth blade 


4^ = 4* 

4> 2 = 4 * + — 
2 2 

4' 3 » 4 ' + 1 1 

%= 4'+— 
2 


The azimuth angle of each blade changes at every instant i.e., it is a function of time 
and is given by, in non-dimensional form, as 
'¥ = Qt= x 


2.4.3 Equations in Non-Dimensional Form: 


The blade equations (2.2) to (2.7) and fuselage equations (2.10) and (2.1 1) are written 
in terms of the real time variable. These equations can be put into non-dimensional 
forms using the non-dimensional time (t ). 

Introducing the following dimensionless quantities 

«=f, T = Ot, =^-0=1,2) 


v(-M) 


v(*>0 

L 


the non-dimensional lag deformation of the blade is assumed to be of the from : 

vg.T) = I (2-13) 

Z=1 


where <(>, is the rotating mode shape of the kth blade, explained in section 2.5. nl is 
the number of rotating modes. For the present study the value of nl is taken as l.The 
blade equations can be re-written as explained below. 


2.4.3.1 Kinetic energy of the blade: 


The kinetic energy of the blade is given in sec. 2.4. LI. Assuming that the blade is 
very stiff in the axial direction, the following inextensionality condition is imposed : 

ii = -If(— ) 2 dx (2.14) 

2' v a* v ' 

On substituting the given inextensionality condition and stated dimenssionless 
quantities, kinetic energy of the blade in non-dimensional from can be written as 
i l x _ 

1 0 0 

+ Q.h6 y cos x ¥ k -Qhdx sin'F*} 2 +{£Q$ 1 (<£)y ljfc (T) + Qx 

Q * - — — 

j{<f>' J (Z)y lit (r)} 2 dx-n/iOysm'i'j t -QhB x cos ^ k } 2 d^ 
2 o 

(2.15) 


2.4.3.2 Strain energy of the blade due to bending: 

The strain energy of the blade is given in sec. 2.4. 1.2, On using eq.(2.3), it can be re- 
written as 

U,=t f(^){f',«)j'i*(r)> 2 <»« ( 2 . 16 ) 

1 o L 

where El is the flexural rigidity of the blade. 


2.4.3.3. Strain energy associated with elastomer: 

The expression for the strain energy of the spring due to deformation A/ is given in 

sec 2.4. 1.3. Using assumed lag mode deformation of the blade it can be written in 
non-dimensional form as 


L /, = I*, r u (r)i 2 -1*3 {I T lt ( T )} 4 +i*5 (i’i'l y,k «> 6 -L 7 {i’i'l n* (V}' 


(2.17) 

Therefore, the spring deformation in lag mode (A/) is written, using assumed lag 
deformation of the blade, as 


A/ *xc{*,« 2 )yi* W)-(«2 - W WH • (2-18) 


or, 



where, 


A^LViYviv) 

*1 4l(h)-($2 -Zl)fl(t2)-il(Z l) (2.19) 

A ^ 

The quantity Ti can be interpreted as normalized deformation of the spring in the 
first lag mode. 

2.4.3 .4. Dissipated energy associated with elastomer: 

In non-dimensional from it can be written as 


W d = {— A/ + Fsgn| A/ 1}.A/ 
(0 

or 


h a'. a' 

W d = {- (OL T, Yu (r)) + F sgn | QL T, Yu (r) |} . {L Ti Yik (t)} 

CO 

Introducing equations (2.15) to (2.20) into variational principle 


( 2 . 20 ) 


5 J(r-!7 t -U,+W d )dt = 0 


(2.21) 


and collecting all the terms, we obtain. 

2 

{-Yik (*) y\k ( T ) - Yu 00 Yu 00} (a,) + {y lk 00 - Yu 00) (a 2 ) + 


(-){0r cos T*y 1Jk (T) -Q x sin T 4 y u (r)}(a 3 ) - {y u (r)}(a 5 ) + 

cos'!-, + 0 , sin'F.Ka.) 




v T a 1 

/-VI/ \4f.,3 


{Tu(t)} + ^-(^) 4 {^(t)} 


mQ 2 f " mQ 2 


jy T9 A / IF T* a * 

W 6 04 (T ) >■ + CP. ) 8 (t ) } 


mQ .03 L 


m 2 M)}+ 


FT, 


mQ 2 L 2 
FT 

^ {?«<*)}(«.) -0 


{sgn|OLT, y,*(t)l} 


mCrL 



where 


i x _ 


X _ 


«1 = J[j»'l«)> 2 *].[/{#'l«)} 2 ]rf? 


0 0 
1 _ 


«2- J»i«)) 2 <*« 

0 

«3 = |»(«)) 2 *]^ 
0 

1 _ 

0 

o, = jf[*]-[ 


0 0 
1 _ 


o« = 


On rearranging the terms, we get four second order coupled non-linear differential 
equations which govern the behavior of the corresponding blade as follows: 

1 


y \k = 


(a 2 +ccij4) 


M){sgn I OL V, y u |} + - (00k **> 


+ (^ ) > + (« 2 - « 5 - a - ^ ) } + (^ ) {^ 3 * > - (y; ) > + (^ ) {^ 7 * > + 

cos T* +0, sin '¥ k }+Ui)$r cos ^ k y lk -6 X sin ^ k y xk )] 

for k=l to 4 (2.22) 

Thus, for the coupled rotor / fuselage model we have a set of six second order non- 
linear differential equations, which have to be integrated in time domain in order to 
analyze the time development of various degrees-of-freedom. 

I 

The non-dimensional coefficients used in eq (2.22) are 

a ' 

F'Pi 
a \ 

mQ. 2 !? 


b x = 


*QPi) 2 

mClcoL 


h =«! 



h\ = 


OCj 


= EI(a 6 ) 
mQ, 2 !? 


rfi- 


e \ 


A ^ 

*ipyi) 2 

mCl 2 L 

a' < 
^iCFiy 
mCl 2 


A 


, a/ , 

g s £ 3 QPi) 6 


mCl* 


/ 


Si = 


/ir 7 z. 3 (y i ) 8 

mO . 2 


m i=(|)(a 4 ) = ^ 


7i =(-)(« 3 ) 

Equations are arranged in the from 


[M] M + [C] M + [AT] M + {Fj^ (y, y, y)} = 0 


(2.23) 


where [A/], [C] and [AT] denote the mass, spring and stiffness matrices, 
respectively. 

The vector {y} represents the degrees of freedom, i.e., 

M = {yn y\2 yn yi4 e x e yf 
Vector { F nl (y, y, y) } includes the non-linear term. 


[ M ] is given by 



-*<*2 0 0 0 ji cos 1 !', p sin MV 

0 -&2 0 0 Pcos^ psinH^ 

0 0 -a 2 0 j3 cos PsinH ' 3 

0 0 0 -a 2 PcosT 4 P sin T 4 

pcos'V l PcosT 2 P cos ¥3 p cos 4^4 -/ t . 0 

P sin 4^ psin^ 2 p sin T 3 p sin T 4 0 -/„ 


Where /3 -Aftm 


(2.24) 


Damping matrix is given by 


— (TO 2 + 

mnwL 


jF'Fi a 4 - * 

-^sgnjQiT, y u | 


m€l l L l 


h A •> 

l^E m + 


mCl 2 ! 2 


■sgn | CILTi y 12 | 


mClwL 

a* 

FT 1 , _ 


(Ti) 2 + 


mn 2 i 2 


{sgniOiT, y 13 | 


0 0 


0 0 


0 0 


2/3siny/| 
- 2/tcos^i 


2 /3 siny/ 2 
- 2/3 cosy/2 


E*3 = 


Stiffness matrix is given by 

a 2 -0C5 -^1 -/1 0 

0 a 2 -a s -d l -li 

0 0 


2/3 sin if/ 1 
-2/3 cos )// 3 


0£ 2 - OC5 - t/j — /j 
0 


fll sgn I niTi y , 4 I 
2/3siny / 4 
- 2/3cosy / 4 

(2.25) 


0 0 


• c, 0 

0 -c„ 


0£ 2 _ « 5 — C?j — /j 0 


(2.26) 
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a 2 - ce 5 -«/, = 




J • J v ' J "’5 „2r -2 _4 mr 1V5/J “>5 

o oo mOrL \mQrU )q 

The vector representing the non-linear term is given by 

lFNtSy*y>y)) = 

f 2 r- 

(«l ^ 11-^11^1) + mo'y cos Vl7l , -6 x sm Wiyi ,} + + (e,)^, 3 ,} - Oi)0^i> + (gi)O'nl 1 

... 2 

■> t a i > M2 y\ 2 - y\ 2 y\ 2 1 + to) { 0 ^ cos y l2 - o x sin 2 } + (Zh {y\2 } + ) 0 & }-(f\){yh)+ (gi ) \y\i } 

- ... 2 1 
(<*1 ) {-yh y\i- yu + <n){Q y «*vsyi 3 -e* sin ^3} + A)^} + (e,)^} - (/i){yf 3 } + (gOiyh) 

2 

t0Ei){-yi4yi 4 ->'i4yi4} + ( T ?){^ cos Wl4-^sinv/4j 1 4} + (^-){3; 1 3 4} + (e 1 ){yi 3 4}-(/i){3 2 i 5 4}+(gl){yi 7 4} 


I Jt _ 


where tj =(-) J[ dx]d$ 


(2.27) 


2.5 Rotating Mode Shape ($, ): 

Using Ref. [5], the comparison function used for the solution is 

£ ^ _ %*[ !LL Mi +3 ) | +3 K 2 , f(i +1 ) ^3j 

1 6 3 6 

Using data given in table (2.2 ) the rotating mode shape is obtained as 

nd A 

/-I 


(2.28) 


(2.29) 


where ND is the number of data points. On normalizing, the rotating mode shape in 
lag motion for the first log mode is given by 

0,(4) = 2.7794294 2 + 2.4865674* -51.6671194 4 + 173.8655404* 

-248.2397714* +119,2092394 7 +36.6896024* + 57.7291324 9 

-275.708214^ 10 +376.621774" -430.8490564 12 + 528.1 136954 * 
-472.7086244 14 + 370.1755024 15 - 670.9403714 16 + 978.5009234 " 
-658.4072244 18 +163.3489794 19 

(2.30) 

This can be used in the blade equations to evaluate the integrals. 
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a l 

-0.054210 

a 2 

- 0.050783 

a 3 

0.372845 

34 

- 0.642042 

*5 

0.315563 

a<s 

0.141461 

a i 

- 0.030847 

as 

-0.216103 

a 9 

0.213805 

aio 

-0.169271 

an 

0.206519 

an 

- 0.197030 

a 13 

0.051855 

a l4 

- 0.126215 

ais 

0.286001 

ais 

- 0.140557 



Table 2.2 Coefficients afs for the rotating mode shape, in 1 st lag mode with 
©„ = 0.7773, from Ref. [6] 
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CHAPTER 3 


SOLUTION PROCEDURE AND TIME 
DOMAIN INTEGRATION 


This chapter describes the methodology used for solving the non-linear ordinary 
coupled differential equations, representing coupled rotor-body dynamics, in time 
domain. The Time Domain Integration approach provides the transient response 
behavior of blades and fuselage with time. 

This chapter is divided into three sections, which deals with, respectively 

i) Solution procedure, for solving the coupled blade and fuselage equations using 
numerical integration technique at every time step. 

ii) Multiblade coordinate transformation. 

iii) Formulation of an index to measure the rate of convergence (stability) or 
divergence (instability) of the system for the prescribed operating and initial 

' conditions, 

3.1 Solution Procedure 

The amplitude-dependent stability analysis of a coupled rotor/fuselage system 
is obtained by simultaneously solving the non-linear coupled equations by numerical 
integration technique in time domain at every time step. 

It can be seen from equations written in the matrix form in section 2.4.3, that 
the mass matrix [M] is time-dependent. Thus, in order to find vector {Y} whi^h 
represents the system degree of freedom the [ M ] has to be inverted at every time 
step. 

The blade and fuselage equations are solved simultaneously by a numerical 
integration routine (NAG)- D02BBF based on Runge Kutta-Marson algorithm for 
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prescribed initial conditions. Blades data are taken from Ref. [1] and are presented in 
table (3.1). 

The response of the system has been evaluated for various combinations of 
rotor rpm and magnitude of initial condition (disturbance). 

To check the validity of the time step, time integration was carried out with 
two time steps viz.. At = 0.0015 sec. and At = 0.003 sec. (real time), both the results 
are found to converge to the same value. Hence, in all subsequent computations, a real 
time step of At = 0.003 sec. is used. The results are produced for £i= 0. 1 and ^ = 0.25. 
In order to study the effect of support damping, the same procedure is repeated one 
with support damping and the other without support damping. 

3.2 Multiblade Coordinate Transformation 

For an N-bladed rotor with blades evenly spaced around, the azimuth angle for the k* 
blade, at any instant, can be written as 

V'A=V'+“(*-0 k- 1,2,3 N ,(3.1) 

where y = Ot, the non-dimensional time variable. 

Let Yik be a generalized coordinate associated with any degree of freedom of 
the k‘ h blade. Since this Yik is associated with the blade, which is rotating, it is called 
a rotating coordinate. If there are N blades then there will be N rotating coordinates. 
By suitably choosing a transformation these N rotating coordinates can be 
transformed into another set of N coordinates, each of which is associated with the 
motion of the rotor disk as viewed from a non-rotating frame. This type of 
transformation is denoted as Multiblade Coordinate Transformation. Basically, this 
transforms the rotating coordinates into a non-rotating frame, which is obtained by the 
following operations. 

1 n : 

Collective operator : Y c = — ^T lt (3*2) 

A' *.j 

Alternating operator : ^ (“!)**»* ( 3 - 3 ) 

A *„i 
1 N 

1-Cosine operator : Y, c ^(cosy k )Y lk (3*4) 

A/ *.i 


25 



I -Sine operator 


I N 

r* =T7 S( sin n)^u 

/V A* l 


(3.5) 


In our case, N = 4. Knowing tiie time response of each blade degree of 
lUvdoin, oik' elm easily ohtniii the time response of tutor degiecs oHieetloiu. 


3.3 Measure Of Instability 

The objective is to formulaic an index, which measures the rate of divergence 
(instability) or convergence (stable response) and to correlate this index with the 
operating and/or initial condition 

Following procedure is adopted, to obtain the index: 

In the divergent case, the non-dimensional time at which the response of the 
blades and/or fuselage becomes double of the given initial disturbance, is taken as the 
index. This is explained in figure [3.1J. 

In the convergent case, the non-dimensional time at which the response of the 
blades and/or fuselage reaches 1/1 0 th of the given initial disturbance, is taken as the 
index. This is explained in figure [3.2]. 

A 3-D plot has been generated with rotor rpm and amplitude as the variables 
in order to study the correlation between the two on the system response. 
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Table 3.1 Blade and non-dimensional fuselage input data for ground resonance 

analysis 
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^napter 4 


RESULTS 


The transient response analysis is carried out to examine the behavior of coupled 
rotor/fuselage system in time, for prescribed operating (rotor rpm) and initial 
conditions. 

The response analysis has been carried for two cases, namely, system with support 
damping and without support damping. The results are presented in the following 
sections. 

4.1 With support damping 

For the baseline configuration the location of elastomer and torque tube are 0.1 
and £, 2 = 0.25, respectively. The damping ratios in the pitch and roll modes are ^q x = 
0.082, Co y =0.0529. 

The results are generated for a range of rotor r.p.m. and initial disturbance in roll 
mode. The r.p.m. range is 250 r.p.m. to 700 r.p.m. The roll input is in the range 
9 X (0)= 0.1 degree to 10 degree. The initial conditions corresponding to other degrees 
of freedom are taken as zero, for all the cases. 

The qualitative nature of response of the system to initial disturbance is shown in 
Table 4.1, The table indicates the region of instability in both r.p.m. and initial 
disturbances in 9 X . The effect of non-linearity of the elastomer is to introduce an 
amplitude dependent instability, which can be clearly seen for the cases corresponding 
to rotor r.p.m., 450-640. At 450 r.p.m., the system is stable for initial disturbance 
0 X (O)< 4 degrees and it is unstable for 9 X (0)> 6 degrees. 

The time response of the individual blade and fuselage degrees of freedom are shown 
in fig. [4.1] to [4.12)] respectively for rotor r.p.m. 250,350,400,450,460,550,640 and 
700. 



Figures [4.1] and [4.2] show the response of all the degrees of freedom in the stable 
zone of operating condition (250,350 r.p.m.) and initial condition, 0 X (O)= 4 degree, 6 
degree, respectively. 

For 450 r.p.m., the results corresponding to four cases 0 X (O)= (0.8deg.,4deg.,6deg. and 
lOdeg.) are shown in figures [4.3] to [4.6]. 

Figures 4.7 to 4.8 shows the response of all the degrees of freedom for 460 rpm and it 
is observed that system is stable for initial disturbance 0 X (O)< 1 degree and it is 
unstable for 0 X (O)> 2 degree. 

I 

Maximum instability is seen for 550 rotor r.p.m., figure [4.9]. 

At 640 r.p.m., the system is stable for initial disturbance 9 X (0)< 4 degree and it is 
stable for 0 X (O)> 6 degree, shown in figure [ 4.10] and [4.1 1]. 

Figure [4.12] shows the response of all the degrees of freedom for operating condition 
700 r.p.m., and initial condition 0 X (O)= 4 degree and it is observed that the system gets 
stable at higher value of initial disturbance also. 


4.2 Without support damping 

For the system without support damping, nature of response of the system to initial 
disturbance is shown in Table 4.2. System response shows a similar behavior except 
for the short range of rotor r.p.m., 350r.p.m. At 350r.p.m., the system is stable for 
initial disturbance 0 X (O)< 4degree and it is unstable for 0 X (O)> 6 degree. 

4.3 Multiblade Coordinate Transfomation 

In . order to study the behavior of rotors degree of freedom, blade degrees of freedom 
has been transformed and the response of blades in collective, alternating cosine- 
cyclic and sine-cyclic modes are presented for following three cases, 

(a) 250 rpm, 0 X (O)= 4 degrees (stable zone) 

(b) 550 rpm, 0 X (O)= 4 degrees (unstable zone) 

(c) 700 rpm, 0 X (O)= 4 degrees (stable zone) 

Figure [4.23]- [4.25] explains the behavior of rotor blades in the transformed 
coordinates. 


4.4 Measure of instability 


Figures [4.1] and [4.2] show the response of all the degrees of freedom in the stable 
zone of operating condition (250,350 r.p.m.) and initial condition, 9 X (0)= 4 degree, 6 
degree, respectively. 

For 450 r.p.m., the results corresponding to four cases 0 X (O)= (0.8deg.,4deg.,6deg. and 
lOdeg.) are shown in figures [4.3] to [4.6]. 

Figures 4.7 to 4.8 shows the response of all the degrees of freedom for 460 rpm and it 
is observed that system is stable for initial disturbance 0 X (O)< 1 degree and it is 
unstable for 9 X (0)> 2 degree. 

Maximum instability is seen for 550 rotor r.p.m., figure [4.9]. 

At 640 r.p.m., the system is stable for initial disturbance 0 X ( 0)£ 4 degree and it is 
stable for 0 X (0)> 6 degree, shown in figure [4.10] and [4.1 1], 

Figure [4.12] shows the response of all the degrees of freedom for operating condition 
700 r.p.m., and initial condition 0 X (O)= 4 degree and it is observed that the system gets 
stable at higher value of initial disturbance also. 


4.2 Without support damping 

For the system without support damping, nature of response of the system to initial 
disturbance is shown in Table 4.2. System response shows a similar behavior except 
for the short range of rotor r.p.m.,350r.p.m. At 350r.p.m., the system is stable for 
initial disturbance 0 X (O)< 4degree and it is unstable for 0 X (O)> 6 degree. 

4.3 Multiblade Coordinate Transfomation 

In order to study the behavior of rotors degree of freedom, blade degrees of freedom 
has been transformed and the response of blades in collective, alternating cosine- 
cyclic and sine-cyclic modes are presented for following three cases, 

(a) 250 rpm, 0 X (O)= 4 degrees (stable zone) 

(b) 550 rpm, 0 X (O)= 4 degrees (unstable zone) 

(c) 700 rpm, 0 X (O)= 4 degrees (stable zone) 

Figure [4.23]- [4.25] explains the behavior of rotor blades in the transformed 
coordinates. 


4.4 Measure of instability 



A 3-D plot has been generated in order to study the correlation between the rate of 
convergence/divergence with the rotor rpm and initial disturbance in roll direction. 
From figure [4.26] it is seen that 

(1) the rate of divergence is maximum around 550 rpm. 

(2) slope of the surface first decreases upto 600 rpm and then starts increasing for 
higher degree. 

A similar plot is shown for stable case also in figure [4.27], 

4.5 Validation 

The result obtained from time domain integration approach matches well with 
the range predicted by the linearized stability analysis. Maximum instability is 
also seen at 550 rpm . A short range of instability occurs at 350 rpm and 
0 X (0)> 6 degree for zero damping case . Results so obtained from the time 
integration approach for the non-linear system matches fairly well with the 
results produced in Ref. [5]. 
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Fig, 4.1(a) Blade response in non-dimensional time domain 
250 rpm, 0x(O) = 4 degree, (^ex = 0.082 , ^qy = 0.0529) 
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Fig. 4.1(b) Fuselage motion in non-dimensional time domain 
250 rpm , 0 X (O) = 4 degree, (Cex = 0.082 , £ 9Y = 0.0529) 
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Fig. 4.2(a) Blade response in non-dimensional time domain 
350 rpm, 0 X (O) » 6 degree, (Cox = 0.082, Coy = 0.0529) 
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Fig. 4.2(b) Fuselage motion in non-dimensional time domain 
350 rpm , 0x(O) = 6 degree, (Cox = 0.082 , Coy = 0.0529) 
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Fig. 4.3(a) Blade motion in non-dimensional time domain 

450 rpm, 9 X (0)= 0.8 degree, (<**=0.082, <*,= 0.059) 
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Fig. 4.3 (b) Fuselage motion in non-dimensional time domain 
450 rpm, 0 X (O)== 0.8 degree, (^ 0x =O.O82, ; Gy =0,0529) 
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Fig, 4.4(a) Blade response in non-dimensional time domain 
450 rpm , 0x(0) “ 4 degree, (Cox “ 0.082 , Coy « 0.0520) 
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Fig. 4.4(b) Fuselage motion in non-dimensional time domain 
450 rpm , 0 X (0) = 4 degree, (fa s 0.082 , fa = 0.0529) 
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Fig. 4.5(a) Blade response in non-dimensional time domain 
450 rpm , 0x(0) = 6 degree, (£qx — 0.082 , £oy = 0.0529) 
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Fig. 4.5(b) Fuselage motion in non-dimensional time domain 
450 rpm , 0x(O) = 6 degree, (^ex = 0.082 , ^qy = 0.0529) 
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Fig. 4.6 (a) Blade motion in non-dimensional time domain 

450 rpm , 0x(O) = 10 degree, (^ex = 0.082 , y = 0.0529) 
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Fig. 4.6 (b) Fuselage motion in non-dimensional time domain 

450 rpm , 0 X (O) = 10 degree, (£ex = 0.082 , ^ eY = 0.0529) 
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Fig.4.6 (c) Motion of blade in non-dimensional time domain 

450 rpm, 0 X (O)= 10 degree, linear elastomer, (<^ x =0.082, £e y =0.0529) 
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Fig.4.6 (d) Motion of the fuselage in non-dimensional time domain 

450 rpm, 6 X (0)= 10 degree, linear elastomer, (£ 0x =O.O82, Cey=0-0529) 
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Fig. 4.7(a) Blade response in non-dimensional time domain 
460 rpm , ©x(0) = 1 degree, (Cex = 0.082 , C,eY - 0.0529) 
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Fig. 4.7 (b) Fuselage motion in non-dimensional time domain 
460 rpm , 0x(O) = 1 degree, (<^ex = 0.082 , ^ev = 0.0529) 
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Kg. 4.8(a) Blade response in non-dimensional time domain 
460 rpm , 8x(0) = 2 degree, (£ex = 0*082 , Cev = 0.0529) 
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Fig. 4.8 (b) Fuselage motion in non-dimensional time domain 
460 rpm , 8 X (0) * 2 degree, (£ex = 0.082 , = 0.0529) 
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Fig. 4.9(b) Fuselage motion in non-dimensional time 
550 rpm , 0x(O) *» 4 degree, (^ox = 0.082 , C,qy = 0.0529) 
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Fig. 4.10(a) Blade response in non-dimensional time domain 

640 rpm , 0x(O) * 4 degree, <£ex = 0.082 , ^ev = 0.0529) 
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Fig. 4.10 (b) Fuselage motion in non-dimensional time domain 

640 rpm , 0x(0) = 4 degree, (£,ex = 0.082 , ^oy = 0.0529) 
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Fig. 4.11(a) Blade response in non-dimensional time domain 
640 rpm , 0x(O) 38 6 degree, (£ex = 0.082 , = 0.0529) 
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Fig. 4.1 1 (b) Fuselage motion in non-dimensional time domain 
640 rpm , Gx(0) = 6 degree, (Cex = 0.082 , ^ev = 0.0529) 
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Fig. 4.12(a) Blade response in non-dimensional time domain 
700 ipm , 0x(O) = 4 degree, (£ex = 0.082 , Cev = 0.0529) 
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Fig. 4.12 (b) Fuselage motion in non-dimensional time domain 
700 rpm , 9 X (0) = 4 degree, (Cex = 0.082 , ^ eY = 0.0529) 
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Fig. 4.13(a) Blade motion in non-dimensional time domain 
250 rpm, 0x(O) = 4 degree, (^ex = 0.0, £ 0Y = 0.0) 
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Fig. 4.1 3(b) Fuselage motion in non-dimensional time domain 
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Fig. 4.14(a) Blade response in non-dimensional time domain 
Operating condition: 350 rpm, Initial condition: 4 degree (£e x =0.0, £e y =0.0) 







8*87 
8*86 
uj 8*85 
§ 8.84 

* 8.83 

d 0* 0 2 

g 0.81 
8 

-8.81 

- 0.02 

8 100 288 388 480 508 600 700 808 900 1908 

HON-DIHENSIONRL TIME 


8.002 

8.0815 

g 8.881 

X 0.0085 

5 0 

g “0.0005 

o- -8.001 

-8.0015 

- 0.002 

0 100 200 300 400 500 600 700 800 900 1000 

NON-DIMENSIONBL TIME 


Fig. 4.14(b) Fuselage motion in non-dimensional time domain 
Operating condition: 350 rpm, Initial condition: 4 degree (£ex = 0.0, Cev = 0.0) 
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Fig. 4.15 (b) Fuselage motion in non-dimensional time domain 

350 rpm, 0x(O) = 6 degree, (^@x = 0.0, C,qy = 0.0) 
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Fig. 4.16 (a) Blade motion in non-dimensional time domain 
450 rpm, 0x ( 0 ) = 4 degree, (£ex = 0 . 0 , C,qy — 0 . 0 ) 
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Fig. 4.16 (b) Fuselage motion in non-dimensional time domain 
450 rpm, 0 X (O) = 4 degree, (£ex = 0.0, £oy = 0.0) 
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Fig. 4.17 (a) Blade motion in non-dimensional time domain 
450 rpm, 0 X (0) = 6 degree, (Cex = 0.0, <^ 0 Y = 0.0) 
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Fig, 4.17 (b) Fuselage motion in non-dimensional time domain 
450 rpm, 9x(0) = 6 degree, (^ox = 0.0, <^oy — 0.0) 
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Fig. 4.18 (a) Blade motion in non-dimensional time domain 
460 rpm, 0 X (O) = 0.8 degree, (^ox = 0.0, 0 >y = 0-0) 
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Fig, 4,18 (b) Fuselage motion in non-dimensional time domain 
460 rpm, 0 X (0) « 0.8 degree, (Cox = 0.0, Coy “ 0.0) 




208 400 600 800 1000 

NOH-DIHENSIONRL TIME 


69 




Y(3? VC2> 


BLADE 1 




BLADE 2 



BLADE 3 



BLADE 4 


Fig. 4.19 (a) Blade motion in non-dimensional time domain 
460 rpm, 0 X (O) = 1 degree, (£gx = 0.0, y = 0.0) 
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Fig. 4.19 (b) Fuselage motion in non-dimensional time domain 
460 rpm, 0 X (O) = 1 degree, (Ce X = 0.0, Coy = 0.0) 
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' Fig. 4.20 (a) Blade motion in non-dimensional time domain 

550 rpm, 0x(0) = 4 degree, (^ox = 0.082, ^oy - 0.0529) 


72 







ROLL MODE 




Fig. 4.20(b) Fuselage motion in non-dimensional time domain 
550 rpm, 0 X (0) = 4 degree, (^ox = 0.0, ^ 0 y = 0.0) 
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Fig. 4.21 (a) Blade motion in non-dimensional time domain 
640 rpm, 0x(O) = 0.8 degree, (£ex = 0.0, y = 0.0) 
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Fig. 4.21 (b) Fuselage motion in non-dimensional time domain 
640 rpm, 0x(O) = 0.8 degree, (<^ox = 0.0, <^oy = 0.0) 
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Fig. 4.22 (a) Blade motion in non-dimensional time domain 
640 rpm, 0 X (O) = 1 degree, (Cox = 0.0, Coy = 0.0) 
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Fig. 4.22 (b) Fuselage motion in non-dimensional time domain 
640 rpm, 0x(O) = 1 degree, (£ex = 0.0, <^qy = 0.0) 
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Fig. 4.23 Multiblade Coordinate Transformation 

250 rpm, 0x(0) = 4 degree, (<^ox — 0.082, <^oy - 0.0529) 
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Fig. 4.24 Multiblade Coordinate Transformation 
550 rpm, G x (0) = 4 degree, (Cox = 0.082, Coy = 0.0529) 
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Fig. 4. 25 Multiblade Coordinate Transformation 

700 rpm, 0x(0) = 4 degree, (i^ox = 0.082, <^oy - 0.0529) 
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Fig. 4.26 Three dimensional plot to show the variation of index of instability with 
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Fig 4.28 Effect of support damping on stability (^ox = 0.082, £oy = 0.0529) 

no support damping (£i= 0.1, %r=- 0.25) 

- with support damping (£i= 0. 1 , 0.25) 
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Chapter 5 


CONCLUSION 


The time domain behavior of the coupled rotor/fuselage system has been studied. The 
properties of the elastomer are taken from the experimental results published in the open 
literature. The major conclusions of the present study are summarised below. 

1. It has been observed that the magnitude of initial disturbance affects the 
aeromechanical stability of the coupled system under ground resonance condition. 

2. The effect of support damping on the response behavior obtained in the present work 
agrees well with that reported in Ref. [6] where lineanzed stability analysis was 
carried out also the range of rotor rpm, in which the system is unstable, matches well 
with range predicted by the linearized stability analysis. 


86 



References 


1 . Huber, H,, “Will Rotor Hubs hose Their Bearings? A Survey of Bearingless 
Main Rotor Development,” Eighteenth European Rotorcraft Forum , Avignon, 
France, Sept, 1992, pp. 506.1-506.19. 

2. Panda, B., and Mychalowycx, E., “Aeroelastic Stability Wind Tunnel Testing 
with Analytical Correlation of the Comanche Bearingless Main Rotor,” 
Journal of the American Helicopter Society, Vol 42, No. 3, 1997, pp. 207-217. 

3. Weller, W. H., “Variation in Hover Aeromechanical Stability Trends with 
Bearingless Main Rotor Design,” Journal of the American Helicopter Society , 
Vol 42, No. 3, 1997, pp. 207-217. 

4. Kampa,K., Enenkl, B., Roth, G. and Polz , G., “ Aeromechanics aspects in the 
design of the EC135” , 23 rd European Rotorcraft Forum, Germany, pp. 38.1- 
38.14, 1997 

5. Pohit, G. Venkatesan, C. and Mallik, A. K., “Influence of Non-linear 
Elastomer on isolated lag Dynamics and rotor/fuselage Aeromechanical 
Stability,” AIAA 2000-1691, 41 st AIAA/ASME/ASCE/AHS/ASC structures, 
structural Dynamics and materials conference, Atlanta, Georgia 

6. Pohit, G. “Dynamics of a Bearingless Helicopter Rotor Blade with a Non- 
linear Elastomeric Constraint, Ph.D Thesis, Dept, of Mechanical Engg., Indian 
Institute of Technology, Kanpur, India, 1999. 

7. Gandhi, F., and Chopra, I., “Analysis of Bearingless Main Rotor Dynamics 
with the Inclusion of and Improved Time Domain Nonlinear Elastomeric 
Damper Model,” Journal of the American Helicopter Society, Vol. 41, No. 3, 
1996, pp. 267-277. 

8. Friedmann, P. P., and Venkatesan, C., “Voupled Helicopter Rotor/Body 
Aeromechanical Stability, comparison of Experimental and Theoretical 
results,” Journal of Aircraft, Vol. 22, No. 2, 1985, pp. 148-155. 


87 



9- Friedmann, P. P., and Venkatesan, C.,” Influence of Unsteady Aerodynamic 

Models on Aeromcchanical Stability in Ground Resonance, " Journal of the 
American Helicoptfe ftolf^l, No. 1, 1986, pp. 65-74. 

10- Bousman, W. G., “An Experimental Investigation of the Effects of Aeroelastic 
Coupling on Aeromechanical Stability of Hingeless Rotor Helicopter,” 
Journal of the American Helicopter Society, Vol. 26, No. 1, 1981, pp. 46-54. 

1 1 - Venkatesan, C., “Influence of Aeroelastic Coupliing on Coupled Rotor/Body 
Dynamics,” 6 th International Workshop on Dynamics, Aeroelastic Stability 
Modelling, of Rotorcraft, University of California, Los Angeles, Nov. 8-10, 
1995. 


88 



